First Measurement of the Fraction of Top Quark Pair Production Through Gluon-Gluon Fusion 
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We present the first measurement of o-{gg tt)/a{pp —* tt). We use 0.96 fb~^ of ^/s =1.96 TeV pp 
collision data recorded with the CDF II detector at Fermilab. We identify the candidate tt events with a high- 
energy charged lepton, a neutrino candidate, and four or more jets. Using charged particles with low transverse 
momentum in events, we find o-(gg —> tt)/a{pp ff)=0.07±0.I4(stat)±0.07(syst), corresponding to a 95% 
confidence level upper limit of 0.33, in agreement with the standard model NLO prediction of 0.15 ± 0.05. 

PACS numbers: 14.65.Ha, 12.38.Qk, 12.38.Aw, 13.85.-t 



Many studies have been dedicated to the understanding of 
the top quark, motivated in part by its large mass that may 
give it a unique role in the generation of mass for the quarks, 
leptons, and force carriers in the standard model (SM) of par- 
ticle physics. In pp collisions at a center-of-momentum en- 
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ergy of y/s =1.96 TeV, (15 ± 5)% of pairs are expected to 
be produced through gluon-gluon fusion and the rest through 
quark-antiquark annihilation based on next-to-leading- 

order (NLO) quantum chromodynamics (QCD) calculations. 
The inclusive ti production cross section has been measured 
by both CDF islHl and DO jsi collaborations using various 
methods and decay modes of the tt pairs, and the results are 
in agreement with SM predictions. However, the details of the 
production process have never been investigated. 

A measurement of the gg tt production cross section 
tests the QCD prediction. Also of interest is any indication 
of new top quark production and decay mechanisms. As the 
partonic cross section calculations are directly related to the 
momentum distributions of constituents of the colliding pro- 
tons such a measurement could assist in reducing the un- 
certainties in the gluon distributions within protons. 

In this Letter we report the first measurement of the frac- 
tional cross section of tt production through gluon-gluon fu- 
sion, (T{gg tt)/a{pp tt). To discriminate between the 
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similar final state signatures of gg ti and qq ti, we 
take advantage of the higher probability for a gluon than for 
a quark to radiate a low-momentum gluon |@]. Therefore, on 
average we expect a larger low-momentum charged particle 
multiplicity in gg — + tt compared to qq — > tt. Given the large 
theoretical uncertainties associated with gluon radiation, we 
do not rely on theoretical calculations for the modeling of the 
charged particle multiplicity. Instead, we use two different 
processes, W^H-n jet and two-jet (dijet) production, with well- 
understood production mechanisms, as calibration samples to 
relate the observed charged particle multiplicity to the fraction 
of processes involving more gluons. 

We use a data sample of y/s =1.96 TeV pp collisions with 
an integrated luminosity of 0.96±0.06 fb^^ recorded by the 
CDF II detector at Fermilab between March 2002 and Febru- 
ary 2006. The CDF II detector is described in detail in ||7t]; 
here, we briefly discuss the components essential for this anal- 
ysis. The detector consists of a tracking system immersed 
in a solenoidal magnetic field of 1.4 T and electromagnetic 
and hadronic calorimeters surrounding the solenoid, followed 
by the muon system. Electrons, photons, and hadronic jets 
are identified using calorimeters and the tracking information. 
Muons are identified by the muon system together with track- 
ing and calorimeter information. The data are collected using 
a three-level trigger system. 

According to the SM top quarks almost always decay to a 
W boson and a bottom quark, and so in t< events we expect to 
have two W bosons and two b quarks. We select tt candidate 
events where one of the W bosons decays to two jets and the 
other decays to a lepton (I) and the corresponding neutrino. In 
this analysis / is either an electron or a muon. Our first calibra- 
tion dataset is a set of li')+njet (n = 0, 1,2, 3) candidate 
events, for which the number of gluons involved in the produc- 
tion process increases with the number of jets |8]. The second 
is a set of events with two back-to-back, high-energy jets. The 
average number of gluons involved in dijet production |9] falls 
with increasing transverse energy (Et) [ 10] of the highest Et 
jet (leading jet), as the relative rate of the qq qq, qg qg 
and gg — + gg subprocesses change. The number of gluons in 
each subprocess is 0, 2, and 4 respectively. 

The W+ n jet data are collected with an inclusive lepton 
trigger that requires an electron with Et >18 GeV or a muon 
with pt > 1 8 GeV/c. We select events with a reconstructed 
isolated electron (muon) candidate with Et > 20 GeV (pT > 
20 GeV/c) and a missing Et (^t)> 20 GeV. We catego- 
rize the W+n jet samples by n, the number of jet candidates 
with Et >15 GeV and pseudorapidity region \ri\ < 2. Jets 



are defined using an iterative cone algorithm 111 lH with a cone 
of Ai? — ^ (A(/))^ + (A77)2 =0.4 and are corrected for abso- 
lute energy response, r\ dependence of calorimeter response, 
and multiple interactions. For the tt data sample, in addition 
to the above, we require four or more jets where at least one 
is identified as originating from a h quark (6 tag). To de- 
fine a h tag, we identify within a jet a long-lived i?-hadron 
candidate through the presence of a displaced secondary ver- 
tex ijsll. In both tt and M^H-n jet samples, we remove any event 



with a second lepton candidate consistent with arising from 
a Z boson decay or a tt event in which both W bosons de- 
cay to leptons. We also veto the events in which the electron 
(muon) is consistent with coming from a conversion photon 
(cosmic ray) I3!]. The dijet data are collected using two inclu- 
sive jet triggers that require a jet with Et of at least 50 GeV 
or at least 100 GeV (Jet50 and JetlOO datasets). We require 
a minimum leading jet Et of 75 and 130 GeV for Jet50 and 
JetlOO datasets, respectively, to avoid any trigger bias. We 
remove events containing an electron (muon) candidate with 
Et >20 GeV (pT >20 GeV/c). We also require exactly two 
jets with Iryl < 2 and a minimum Et of 20 GeV and with the 
two jets back-to-back, having a |A0| > 2.53 rad. 

The background processes in our tt sample consist of W+ 
jets, electroweak processes (WW, WZ, ZZ), single top 
quark, and multi-jet QCD processes (non-W). For non-M^ 
and W+ jets background, we can have a real b tag (heavy 
flavor background, HF) or have a b tag due to misidentifica- 
tion (light flavor background, LF). We estimate LF and HF in 
events with a real W boson using various calibration datasets. 
For the small fraction of events from non-VK sources, we as- 
sume the non-T4^ background is equal parts HF and LF. The 
results of the analysis are insensitive to this assumption. Sin- 
gle top quark processes are part of HF, while diboson back- 
grounds, ignoring the few Z ~^ bb events, are included in LF. 
We find 240 tt candidates with an estimated background con- 
tamination of (13±2)%. The background estimates are found 
using the method explained in I3j] . 

The number of low-pT charged particles A'trk is affected by 
low-energy particles arising from jet fragmentation as well as 
multiple interactions within the same pp bunch crossing. To 
include a track in our definition of A^trk, we require it to have 
apT in the range 0.3-2.9 GeV/c and \r]\ < 1.1 and to originate 
from the vertex associated with the charged leptons and jets. 
We reject the track if it fafls within AR =0.6 and AR =0.4 of 
jets with Et >15 GeV (high Et jets) and 6< Et <15 GeV 
(low Et jets), respectively. Excluding these tracks results in 
different available tracking area for each event. We therefore 
correct the observed multiplicity to the total tracking coverage 
in 77 and cj) event-by-event. The resulting track multiplicity still 
has a modest dependence on the number of high Et jets in 
the event. We therefore make a further correction to iVtrk by 
measuring this dependence in multi-jet QCD candidate events 
and using this as a per-jet correction(^ 1 track per jet) to the 
multiplicity for all jets with |?7| < 1.1. 

We show that there is a correlation between the average 
number of low-pT charged particles (A^tik) and the average 
number of gluons involved in the production process (Ng) in 
a given sample. We count the number of gluons that are part 
of the production process using Monte Carlo (MC) calcula- 
tions for both dijet and the W+ n jet data samples. The W+ n 
jet MC sample is created using the ALPGEN 1 12] program fol- 
lowed by PYTHIA L13] to perform the jet fragmentation. The 
MC dijet events are created using the PYTHIA MC. We plot 
the observed (iVtrk) in data against the expected (Ng) from 
MC calculations for the calibration samples in Fig. [T] This 
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FIG. 1: The correlation between the average number of 
low-pT charged particles (data) and the average number of 
gluons (MC). The dotted line is from a linear fit to the points. 



FIG. 2: The number of low-px charged particles for the tt 
candidates, the fit result, the gluon-rich, and no-gluon 
components. 



demonstrates an approximately linear dependence between 
(-^trk) and (Ng). We do not use this plot to obtain our re- 
sult, but rather directly fit the observed iVtrk distributions as 
described below. 

The (A^trk) and (Ng) correlation enables us to define A'trk 
distributions each representing a specific average number of 
gluons involved in the production process. We use this corre- 
lation and the observed TVtrk distributions in the W+ jet sam- 
ple and the dijet sample with leading jet Et of 80-100 GeV 
to define a no-gluon and a gluon-rich A'trk distribution, re- 
spectively. To do so, an iterative procedure is adopted in or- 
der to remove the no-gluon (gluon-rich) contribution from the 
iVtik distribution of the 80-100 GeV dijet (W+ jet) sam- 
ple. There are no significant changes in the distributions after 
the first iteration. The W+ jet sample is largely composed 
of the Drell-Yan qq' process with a small QCD background 
of order 4% and contribution from W production in associa- 
tion with other partons where none of the final state jets are 
detected. The fraction of W+ jet candidates with produc- 
tion processes involving gluons is estimated to be (5 ± 4)%. 
The no-gluon contribution of dijet candidates with leading jet 
Ex of 80-100 GeV comes from qq qq processes and is 
estimated to be (27 ± 3)%. 

To verify that the no-gluon or gluon-rich distribution can 
model the A'trk distribution of any process with comparable 
(Ng) regardless of the center-of-momentum energy, we check 
the W+ 1 jet data sample, and we see no dependence on jet 
Erin (iVtrk). 

The gluon-rich fraction associated with a given iVtrk distri- 
bution can be found using a binned UkeUhood fit of the A'trk 
distribution of the form 



NifgJ'giNtrk) + (1 - fg)^q{Ntrk)l (D 



where N is the normalization factor and one of the free pa- 
rameters, fg is the fraction of gluon-rich components of the 
sample and the second free parameter, and .^-^(A'trk) and 
jFg(7Vtik) are the normalized gluon-rich and no-gluon param- 
eterizations, respectively. We have verified this technique is 
free of bias using randomly generated samples from the data 
iVtik distributions. 

The iVtrk distribution of the candidates, shown in Fig.|2] 
has a mean of 10.6±0.5. The fit, shown in the figure, mod- 
els the data distribution very well, based on a goodness of 
fit test with 92% probability. The measured gluon-rich frac- 
tion in tt candidates determined by fitting the A'trk distribu- 
tion consists of two components, the tt gluon-rich fraction 
and the background gluon-rich fraction. Therefore, knowing 
the background fraction in our sample and the measured fg 
from the fit, we can write 



fg = fbC + (1 - fb)ff, (2) 

where f^^^ and /** are the gluon-rich fraction of the back- 
ground and signal, respectively. To estimate fg^^, we mea- 
sure fg in the W+ 1, W+ 2, and W+ 3 jet data samples with 
no b tag and with at least one b tag using the fit to the Nf^k 
distribution for each of these samples. We then extrapolate 
the fg values from the W+ 1, 2, and 3 jet samples to W+ 4 
or more jet bins for both &-tag and no &-tag samples. We con- 
sider the 6-tag sample as representative of HF and the no b- 
tag sample as representative of LF. Using these extrapolations 
and the LF and HF fractions, we find f^^i^ = 0.54 ± 0.09 and 
/f = 0.09 ±0.16. 

Given /**, we measure (T{gg — > tt)/a{pp tt) as 
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= 0.07±0.14(stat), 



where Agg and Aqq are the acceptance for gg 



_ (3) 
tt and 

qq tt, respectively. Using PYTHIA [13fl MC calculations, 
we find (14.1 ± 0.5)% and (11.5 ± 0.4)% for Agg and Aqq, 
respectively. The acceptance uncertainties include the system- 
atic uncertainties. 

The systematic uncertainties of this measurement, a total of 
0.07, arise from uncertainties in the measurement of iVtrk and 
the subsequent calculations. The uncertainties in A^trk are due 
to the per-jet correction(0.05), the estimated gluon content of 
the W+ jet sample(0.04), and the choice of the low Et jet 
cut(0.02). In addition to these sources, there are uncertainties 
associated with the estimated qq qq fraction of the 80- 
100 GeV dijet sample, the background fraction, the modeling 
of the background gluon-rich fraction, the non- background 
fraction, and the acceptances; these are all negligible. To esti- 
mate the effects of all the above uncertainties, we changed the 
central values and measured the change in the relevant vari- 
ables. 

The result corresponds to an upper limit of 0.33 at 95% 
confidence level. We use a classical statistical technique to set 
the limit by simulating the possible outcomes for a given true 
value taking into account the systematic effects. 

In conclusion, we have presented the first measure- 
ment of (j{gg tt)/a{pp tt) and found 
0.07±0.14(stat)±0.07(syst), corresponding to an upper limit 
of 0.33 at 95% confidence level, in 0.96 fb^^ of data col- 
lected at CDF. This is in agreement with the SM prediction 
of 0.15 ± 0.05, and does not suggest that non-standard model 
processes 11411 contribute to top quark pair production at the 
Tevatron. 
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